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ABSTRACT 

We present new deep optical spectra of 9 high-z radio galaxies (HzRGs) at z > 2.7 obtained with F0RS2 on VUT. These rest-frame 
ultraviolet spectra are used to infer the metallicity of the narrow-line regions (NLRs) in order to investigate the chemical evolution 
of galaxies in high-z universe. We focus mainly on the C iv/He ii and C iii]/C iv flux ratios that are sensitive to gas metallicity and 
ionization parameter Although the N v emission has been widely used to infer the gas metallicity, it is often too weak to be measured 
accurately for NLRs. By combining our new spectra with data from the literature, we examine the possible redshift evolution of the 
NLR metallicity for 57 HzRGs at 1 S z S 4. Based on the comparison between the observed emission-line flux ratios and the results 
of our photoionization model calculations, we find no significant metallicity evolution in NLRs of HzRGs, up to z ~ 4. Our results 
imply that massive galaxies had almost completed their chemical evolution at much higher redshift (z > 5). Finally, although we 
detect strong N v emission lines in 5 HzRGs at z > 2.7, we point out that high N v/He ii ratios are not indicative of high metallicities 
but correspond to high ionization parameters of gas clouds in NLRs. 

Key words, galaxies: active - galaxies: evolution - galaxies: nuclei - quasars: emission lines - quasars: general 
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1. Introduction 

Chemical evolution of galaxies is one of the most important 
aspects to understand the formation and evolution of galaxies, 
since it is closely related with the past star formation history of 
galaxies. The most straightforward way to investigate the chem- 
ical evolution of galaxies is measuring the metallicity of galax- 
ies at various redshifts and exploring the systematic trends of the 
metallicity as a function of redshift (see Maiolino et al. 2008 and 
references therein). The gas-phase metallicity of star-forming 
galaxies in the local universe can be measured rather easily by 
analyzing optical emission line spectra (e.g., [On]/i3727, H/3, 
[Om]/15007, Ha, [Nii]^6584; see Nagao et al. 2006b and refer- 
ences therein). However, this is difficult at z > 1, since these 
emission lines are generally very faint and shifted into near- 
infrared wavelengths where sensitive spectroscopic observations 
are more difficult than in the optical. 

An alternative approach is to focus on active galactic nu- 
clei (AGNs) instead of star-forming galaxies. AGNs generally 
show various emission lines in rest-frame ultraviolet to infrared 
wavelengths, that arise in gas clouds photoionized by the ra- 
diation from their central engines. Here we focus on the rest- 
frame ultraviolet lines (e.g., Nv/11240, Civ /11549, Heii/11640, 
Cm] /1 1909), since we can easily measure the emission-line 
fluxes of high-z AGNs by means of optical spectroscopic obser- 
vations. These UV emission lines in AGN spectra are generally 



* Based on data obtained at the VLT through the ESO programs 
076.A-0664 and 077.B-0051 



brighter than optical emission lines in star-forming galaxies at 
similar redshifts. 

The metallicity of AGNs has been studied extensively, espe- 
cially by focusing on the broad-line regions (BLRs). It has been 
reported that the BLR metallicity (Zblr) is typically higher than 
the solar value (e.g., Hamann & Ferland 1992; Dietrich et al. 
2003; Nagao et al. 2006c), reaching as much as Zblr ~ ISZ© 
in extreme cases (Baldwin et al. 2003b; Bentz et al. 2004). 
However it is not clear how the gas metallicity inferred from the 
broad lines is related to the chemical properties of the host galax- 
ies, since broad lines of AGNs sample only very small regions in 
galactic nuclei (Rblr < 1 pc; e.g., Suganuma et al. 2006), which 
may have evolved more quickly ffian their host galaxies. For this 
reason, here we focus on narrow-line regions (NLRs), instead of 
BLRs. 

Contrary to BLRs, the size of NLRs is roughly comparable 
to the size of their host galaxies (/?nlr ~ 10^"^ pc; e.g., Bennert 
et al. 2006). The mass of the NLR is about Mnlr ~ 10^"^ Mq 
which is much larger than that of the BLR (Mblr ~ 10^-10"* 
Mq, see, e.g., Baldwin et al. 2003a) and the velocity dispersion 
of the NLR emission lines well traces the kinematics of their 
host galaxies (e.g., Greene & Ho 2005; Bennert et al. 2006). The 
NLR is therefore a good tracer of chemical properties on the 
galactic scales. 

In order to study NLR metallicities it is necessary to focus 
on type-2 AGNs, where BLR and strong continuum emission are 
obscured and do not alfect the accuracy of narrow line emission 
measurements. However only very few optically-selected high- 
z type-2 QSOs have been discovered so far. We then focus on 
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Table 1. Journal of Observations 



Source 




E{B - 


Exp." 


Date 


TN J0121+1320 


3.516 


0.039 


23400 


2005 Oct 28, Nov 20, 










2006 Sep 17, 20 


TN J0205+2242 


3.506 


0.100 


5400 


2005 Oct 1 


MRC 0316-257 


3.130 


0.014 


10800 


2006 Jul 27, Aug 25, Oct 2 


USS 0417-181 


2.773 


0.048 


2700 


2006 Oct 2 


TN J0920-0712 


2.760 


0.041 


10800 


2006 Apr 3, 4 


WN J1123+3141 


3.217 


0.018 


14400 


2005 Dec 28, 










2006 Apr 5, May 27, Jun 20 


4C 24.28 


2.879 


0.018 


10800 


2006 Apr 23 


USS 1545-234 


2.755 


0.257 


14400 


2006 Apr 5, 23,24 


USS 2202+128 


2.706 


0.067 


10800 


2006 Jul 2, Aug 16 



" The redshift based on NED (NASA/IPAC Extragalactic Database). 
* Galactic extinction [mag] given by Schlegel et al. (1998). 
" Total exposure time [sec]. 



NLRs of high-z radio galaxies (HzRGs) since several of them 
have been identified even at high-z. 

Some studies on the metallicity of the NLR in HzRGs 
have already been carried out in the past. By studying 
the emission-Une flux ratios of Nv/I1240/Crv/I1549 and 
Nv/I1240/He 11/11 640, generally used to measure the metallic- 
ity of BLRs in AGNs, De Breuck et al. (2000) found that the 
typical metallicity of HzRGs is roughly O.4Z0 < Znlr < S.OZ©. 
They also reported a possible metaUicity evolution in their sam- 
ple: radio galaxies at z > 3 have lower metallicity than those at 
z < 3. However, the emission-line flux of N v in narrow-line ra- 
dio galaxies is generally too faint, especially for metal-poor gas. 
Therefore, only upper-limits on N v are available for the major- 
ity of the radio galaxies investigated by De Breuck et al. (2000). 
Iwamuro et al. (2003) found lower metallicities in the NLRs in 
HzRGs at 2.0 < z < 2.6 than De Breuck et al. (2000), based 
on their analysis of rest-frame optical emission lines from their 
sensitive near-infrared observations. In contrast, Humphrey et al. 
(2008) recently reported that there is no significant difference in 
NLR metallicities between their sample of HzRGs at z ~ 2.5 and 
the lower-zradio galaxies investigated by Robinson et al. (1987). 
To solve these apparent contradictions, it is mandatory to carry 
out further spectroscopic observations of HzRGs. 

Nagao et al. (2006a) proposed a new metallicity diagnos- 
tic diagram that consists of the CiV/11549, Hen/11640 and 
Cm] /1 1909 emission lines, all of which are moderately strong 
in rest-frame ultraviolet spectra of HzRGs, even at low metal- 
hcities. They studied the metalhcity of NLRs of HzRGs and re- 
ported that the observational data do not show any evidence for 
a significant evolution of the gas metallicity in NLRs within the 
redshift range 1 .2 < z < 3.8, and instead they found a clear trend 
for more luminous AGNs to have more metal-rich gas clouds. It 
should be noted, however, that their sample includes only 5 ob- 
jects at z > 2.7 (or only 2 objects at z > 3.0). Thus, observing 
more HzRGs is crucial to assess the possible metaUicity evolu- 
tion. 

In this paper we report new spectroscopic observations of 
9 HzRGs at 2.7 < z < 3.5. By combining the new data with 
the Nagao et al. (2006a) database, we discuss the chemical evo- 
lution and the metallicity-luminosity relation of HzRGs in the 
z ~ 1 - 4 redshift range. We adopt a concordance cosmology 
with (Dm, 0.a) = (0.3, 0.7) and Ho = 70 km s ' Mpc'^ 



2. Observations 

We observed 9 HzRGs at z > 2.7 with FORS2 (FOcal Reducer 
and low dispersion Spectrograph 2) at the VLT (Very Large 
Telescope). These targets were selected from the catalog of De 
Breuck et al. (2000), excluding HzRGs whose C rv. He n and 
C III] have been already measured. All the observations were ex- 
ecuted in the service mode. The list of the target objects with the 
observation log is given in Table 1 . 

Observations were performed with the 300V dispersion el- 
ement and the G435 order-sorting filter to cover the 4450 A < 
/lobs < 8700 A range, which, at the redshift of our targets, cor- 
responds to the rest-frame ultraviolet and thus includes the UV 
emission lines needed for abundance determination. We adopted 
an on-chip binning mode of 2 x 2 with a spectral dispersion of » 
3.2 A pixeL^ and a spatial scale along the slit of » 0'.'25 pixel"' . 
The slit width was I'.'O, narrower than the typical seeing (FWHM 
~ r.'4). Individual exposure time was 900 s, and the total inte- 
gration time for each target is given in Table 1 . 

Standard data reduction procedures were performed used 
available IRAF tasks. Bias was subtracted by using an averaged 
bias image and flat-fielding was performed using dome-flat im- 
ages. Cosmic -ray events were then removed. The wavelength 
calibration was done by using sky lines. After sky subtraction, 
we extracted spectra using apertures of 2'.'25 along the slit. The 
spectra were then flux-calibrated by using the following spec- 
trophotometric standard stars. Note that the weather condition 
was good during most of our observations. Examination of the 
ESO sky monitor database indicates that most (> 90 %) of our 
exposures were taken under photometric conditions. Small vari- 
ations in the absolute calibration due to varying seeing and trans- 
parency are unimportant given that we bin our data into three 
broad luminosity bins. 

3. Observational results 

The final reduced spectra of the observed HzRGs are shown in 
Figure 1. In all spectra continuum emission is very faint; typ- 
ically S/N~ 3 and ranging 1 ^ S/N ^ 10 at Arest ~ 1700 A. 
However the equivalent widths of emission fines are very high 
and their S/N is accordingly very large. The composite spectrum 
of HzRGs, which is shown in the lower right panel of Figure 1, 
was obtained by stacking the 9 spectra of HzRGs after convert- 
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Fig. 1. Spectra for 9 objects observed with FORS2 and the composite spectrum of the HzRGs (lower right), improving by biiming 
x2 in wavelength direction. Each spectrum is shown on two different flux scales, the first to show strong emission Unes and the 
second to show continuum emission and weak emission lines. 
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Table 2. Emission-Line Measurements 



Source 


z" 


Line 


^obs 


Flux* 


FWHM" 








(A) 


(10""^ erg cm-2) 


(A) 


TN J0121 + 1320 


3.517 


Nv^l240 


— 


< 0.094'' 


— 






CiV/11549 


6996.9 


0.263 ± 0.005 


34.4 






Heii/il640 


7414.3 


0.330 ± 0.012 


20.1 






Cm] /1 1909 


8622.2 


0.282 ± 0.009 


32.6 


TN J0205+2242 


3.507 


Nv/il240 


— 


< 0.097"* 


— 






CiV/11549 


6981.8 


0.873 ± 0.025 


33.0 






Heii/il640 


7387.5 


0.519 ± 0.046 


47.3 






Cm] /1 1909 


8586.0 


0.418 ± 0.049 


33.4 


MRC 0316-257 


3.130 


Nvil240 


— 


< 0.038'' 


— 






CiV/11549 


6396.8 


0.267 ± 0.011 


52.4 






He 11^ 1640 


6771.5 


0.301 ± 0.009 


27.9 






Cm] /1 1909 


7874.1 


0.345 ± 0.018 


41.1 


USS 0417-181 


2.773 


Nv^l240 


- 


< 0.103"* 


- 






CiV/11549 


5845.1 


0.356 ± 0.026 


37.5 






He 11^1640 


6188.5 


0.492 ± 0.019 


28.9 






Cm] /1 1909 


7192.8 


0.553 ± 0.047 


30.7 


TN J0920-0712 


2.758 


Nv^l240 


4663.9 


1.015 ± 0.014 


31.6 






CiV/11549 


5821.4 


3.365 ± 0.010 


26.6 






He nil 640 


6161.5 


2.063 ± 0.011 


21.4 






Cm] /1 1909 


7163.2 


1.945 ± 0.028 


28.3 


WN J1123+3141 


3.221 


Nvil240 


5238.7 


1.698 ± 0.013 


31.7 






CiV/11549 


6538.8 


1.570 ± 0.011 


24.8 






He nil 640 


6916.9 


0.425 ± 0.014 


19.9 






Cm] /1 1909 


8046.6 


0.183 ± 0.028 


20.3 


4C 24.28 


2.913 


Nvil240 


4853.2 


1.225 ± 0.012 


46.4 






CiV/11549 


6061.6 


1.235 ± 0.020 


50.7 






He nil 640 


6416.3 


0.978 ± 0.011 


29.5 






Cm]/11909 


7463.8 


0.812 ± 0.041 


41.9 


USS 1545-234 


2.751 


Nvil240 


4642.9 


1.335 ± 0.031 


50.1 






Civil 549 


5810.3 


1.343 ± 0.021 


24.7 






He Hi 1640 


6150.7 


0.878 ± 0.012 


16.9 






Cm]il909 


7150.9 


0.606 ±0.031 


26.9 


USS 2202+128 


2.705 


Nvil240 


4595.0 


0.160 ±0.019 


30.4 






Civil 549 


5739.2 


0.704 ± 0.012 


22.9 






He Hi 1640 


6074.3 


0.289 ± 0.010 


20.3 






C III] i 1909 


7065.5 


0.292 ±0.011 


25.6 



" Redshift of the target calculated from the observed C iv wavelength. 

Line fluxes corrected for the Galactic reddening. 
^ Observed FWHM (includes instrumental broadening. 

3cr upper-limit flux. 



ing to rest frame, normalizing to the He ii flux, and performing a 
3 sigma-clipping rejection. 

Fluxes of observed emission lines are given in Table 2, where 
the measurements were performed by fitting a Gaussian profile 
with the IRAF task splot. The flux error was estimated based 
on the RMS of the pixels around each emission line. The de- 
tected emission lines are well spectrally resolved and instrumen- 
tal broadening is negligible in most cases. Observed emission- 
line widths without correction for instrumental broadening are 
also presented in Table 2. We also provide for each HzRG the 
redshift calculated from the C iv line. For the composite spec- 
trum the flux ratios of Civ/Heii, Ciii]/Civ, and Nv/Civ are 
1.618, 0.628, and 0.584, respectively. 

All the emission lines of C iv. He ii and C iii] (required for the 
metallicity diagnostic method proposed by Nagao et al. 2006a) 
are detected in the spectra of all 9 HzRGs observed. In addi- 
tion, a significant Nv emission is also detected in 5 out of 9 
objects. We provide a 3cr upper limit on the N v emission for ob- 
jects without a significant detection, by assuming an emission- 



line width of ~ 10 A in rest frame. As shown in Figure 1, vari- 
ous fainter emission lines are also detected including Si n/11265, 
Oi-HSi 11^1305, C 11^1335, Oiv]/I1402, Oiii]^1663, Siii/11808, 
[Ne iv]/l2422. We will discuss the properties of these faint emis- 
sion lines in a forthcoming paper (Matsuoka et al. in prep.). 

4. Photoionization model 

To infer the metallicity from the observed emission-line spectra, 
we carried out model calculations by using the photoionization 
code Cloudy version 07.020 (Ferland et al. 1998). Here we as- 
sume that the clouds in the NLR of HzRGs are mainly photoion- 
ized and not significantly affected by shocks. Although Nagao 
et al. (2006a) demonstrated that this assumption is appropriate 
when focusing on C iv. He ii and C iii], we will examine how this 



' We proved for some models that the results of the calculations do 
not change significantly with the use of Cloudy version 08.00 instead 
of version 07.02; the dilference is ~ 10 % at most. 
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Fig. 2. Predicted emission-line luminosity normalized by the 
number of H-ionizing photons, as functions of the NLR metallic- 
ity. The results for N v, C iv, He ii, and C iii] are shown, from top 
to bottom. Red and Green lines denote the models with hydro- 
gen density of hh = lO^cm"^ and lO^cm"^, respectively. Solid 
and dashed lines denote the models with ionization parameter of 
U - 10"'-'' and 10"^ °, respectively. 



assumption is valid in section 5.1. The parameters for the cal- 
culations are (1) the spectral energy distribution (SED) of the 
photoionizing continuum radiation; (2) the hydrogen density of 
a cloud (hh); (3) the ionization parameter (U), i.e., the ratio of 
the ionizing photon density to the hydrogen density at the irradi- 
ated surface of a cloud; (4) the column density of a cloud (A^h); 
and (5) the elemental composition of the gas. 

For the SED of the ionizing photons, we used the "table 
AGN" command, that roughly reproduces the typical SED of 
ionizing photons in AGNs (Mathews & Ferland 1987). We do 
not examine the dependences of the calculations on the SED but 
adopt only this SED, because Nagao et al. (2006a) already in- 
vestigated the dependences in HzRGs and concluded that the 
SED effects do not affect the discussion of the NLR metallic- 
ity evolution in HzRGs. We adopted gas clouds with the hy- 
drogen density rtu = 10^ and 10"* cm"^ and the ionization pa- 
rameters U - lO"^"* - 10"''^. Here we assumed dust-free gas 



clouds, since dusty models are inconsistent with observations 
when high-ionization emission lines of HzRGs are concerned 
(Nagao et al. 2006a). Note that this is suggested also by rest- 
frame optical or near-infrared emission lines (e.g., Marconi et al. 
1994; Ferguson et al. 1997b; Nagao et al. 2003). For the chemi- 
cal composition of gas, we assumed that the all metals scale by 
keeping solar ratios except for He and N. For helium, we as- 
sumed a primary nucleosynthesis component in addition to the 
primordial value. Nitrogen scales as the square power of other 
metal abundances because it is a secondary element. We adopted 
the analytical expressions for the helium and nitrogen relative 
abundances as functions of the metallicity given in Dopita et al. 
(2000). Another free parameter in our calculation is the cloud 
column density. Since we are now focusing on relatively high- 
ionization emission lines (N v, C iv. He n, and C iii]), we stop our 
calculations when the hydrogen ionization fraction drops below 
15%. This requirements ensures that the line fluxes of interest 
will not depend on the choice of a particular column density. 

In Figure 2 we show the dependence of the emission-line lu- 
minosity on the NLR metallicity, for models with nn = 10^ cm"^ 
and 10^ cm and U = 10"'^ and 10"^". The line luminosity 
is normalized by the number of H-ionizing photons in the input 
continuum emission. The model behavior is completely different 
between He ii and the other emission lines: the He ii is a recombi- 
nation line and its luminosity is proportional to the rate of He^- 
ionizing photons, while the other emission lines are collisionally 
excited and their emissivity strongly depends on the gas tem- 
perature. The equilibrium temperature of ionized gas clouds is 
sensitive to metallicity because of the radiative cooling by metal 
emission lines. Therefore the luminosity of collisionally excited 
emission lines decreases at high metallicity. Note that the Nv 
luminosity decreases at much higher metallicity than the Civ 
and C III] luminosity. This is because the nitrogen abundance in- 
creases at high metallicity, i.e., N/H oc (O/H)^. All of the above 
results are insensitive to the adopted gas density and ionization 
parameter. 

Figure 2 suggests that the Civ/Hen flux ratio is possibly 
good metallicity diagnostic. However, this flux ratio also de- 
pends on other parameters such as the ionization parameter. This 
degeneracy can be solved by combining it with the C iii]/C iv 
ratio which is primarily sensitive to the ionization parameter. 
Therefore a diagram involving both the C iv/He ii and C iii]/C iv 
ratios is expected to be a powerful metallicity diagnostic, as orig- 
inally proposed by Nagao et al. (2006a). In Figure 3, the re- 
sults of our photoionization model calculations are plotted on 
this diagnostic diagram. The model grids indicate that this di- 
agram is quite useful to investigate the NLR metallicity. Note 
that, as shown in Figure 3, the C iv/He ii flux ratio depends also 
on the gas density and thus the gas metallicity is not uniquely 
determined through this diagram. However, it can be useful to 
investigate relative differences in metallicity and to assess the 
possible redshift evolution of the NLR metallicity. In the rest of 
the paper we only consider models with riu = lO'^cm"-' (a typical 
NLR density; see, e.g., Nagao et al. 2001a, 2002a). 



5. Discussion 

5.1. The metallicity-redshift and metallicity-lurninosity 
relations in HzRGs 

By comparing the observed flux ratios of HzRGs with the predic- 
tion of the photoionization models, we investigate the gas metal- 
licity. In Figure 4, we plot the flux ratios of HzRGs on the diag- 
nostic diagram with the calculated model grids. We also plot the 
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Fig. 3. Calculated model grids plotted on a diagram of Civ/Hen versus Cm]/Civ. Models for gas clouds with nn = 10^ and 
10'* cm~^ are presented in the left and right panels, respectively. Constant metaUicity and constant ionization parameter sequences 
are denoted by soUd and dashed lines, respectively. 



flux ratios of the composite spectrum of HzRGs on this diagram. 
Though new 9 data are at higher redshift than the data in Nagao 
et al. (2006a) on average, there is no systematic difference in the 
data distribution between our new observations and the data of 
Nagao et al. (2006a). This result naively suggests that there is 
no significant chemical evolution in NLRs of HzRGs, even up to 
z~4. 

We now investigate whether shock models can explain the 
observed flux ratios of HzRGs. In Figure 4 we overplot the shock 
models and shock plus precursor models presented by Allen et 
al. (2008). Specifically, we plot shock models with shock veloc- 
ities of 200 km s"' < Vj < 900 km s"', magnetic parameters 
of Bjn^l'^ = 0, 1, 2, and 4 jxG cm^^^, and solar abundance, and 
shock plus precursor models with shock velocities of 200 km 
s"' < Vj < 400 km s"', a magnetic parameter of B/n'^^ = 2 jjG 
cm-'^^, solar abundance. As shown in Figure 4, both pure shock 
models and shock plus precursor models can not explain the ob- 
served distribution of line ratios, at variance with pure photoion- 
ization models. We therefore conclude that the dominant ioniza- 
tion mechanism of gas clouds emitting C rv. He n, and C ra] is 
photoionization. 

To investigate the relationship between metallicity and red- 
shift or luminosity, we divided our sample into the following 
three redshift bins: 1.2 < z < 2.0, 2.0 < z < 2.5, and 
2.5 < z < 3.8. As for the luminosity, since the non-thermal 
AGN continuum is obscured in HzRGs, we used the He n line 
luminosity as a proxy for the AGN's bolometric luminosity. This 
directly traces the rate of He'^-ionizing photons which is propor- 
tional to the bolometric luminosity. Luminosity bins are as fol- 
lows: 41.5 < logL(Hen) < 42.5, 42.5 < logL(Hen) < 43.0 and, 
43.0 < log L(He n) < 45.0 (see. Table 3). In Figures 5 and 6, we 
investigate the averaged flux ratios of C iv/He ii and C m]/C iv for 
each luminosity and redshift bin to examine the possible metal- 
licity dependence on luminosity and redshift. There is no clear 
dependence between the flux ratios and redshift as shown in 
Figure 5, suggesting that the NLRs metallicity of HzRGs does 
not show significant redshift evolution in the range of 1 5 z ^ 4. 



On the other hand, in Figure 6 we find that the flux ratios are 
correlated with the luminosity, suggesting that more luminous 
HzRGs have more metal-rich gas clouds. Figure 7 shows the 
redshift evolution of average metallicity as derived from our 
sample. Here we use only the HzRG samples with logL(Heii) 
= 42.5 - 43 .0 to minimize possible Malmquist bias. Note that the 
metallicity shown here are derived by assuming nn = 10"* cm"^ 
(see Table 3). Figure 7 clearly shows that there is no significant 
decline in the NLR metallicity at high redshift. In the same way. 
Figure 8 shows the luminosity dependence of the averaged NLR 
metalhcity in our sample. Figure 8 clearly shows the correlation 
between the AGN luminosity and the NLR metaUicity. These re- 
sults are completely consistent with the results of Nagao et al. 
(2006a) but more reliable since our new sample includes a larger 
number of HzRGs at high redshift (z > 2.7). Note that the posi- 
tive correlation between the AGN luminosity and the gas metal- 
licity of ionized clouds is indicated also for the BLR, that does 
not show any redshift evolution (e.g., Hamann & Ferland 1993; 
Nagao et al. 2006c; see also Jiang et al. 2007; Juarez et al. 2009). 

5.2. Interpretation 

In section 5.1, we reported a significant positive correlation be- 
tween the NLR metalhcity and the AGN luminosity. In this sec- 
tion we discuss two possible origins for this correlation: 

- the galaxy mass-metaUicity relation (assuming Lagn °^ 

Mbh oc Mhost) 

- a dependence of NLR metaUicity on the Eddington ratio 

We discuss each of these scenarios below. 

One possible scenario for the positive correlation between 
NLR metallicity and AGN luminosity is that it reflects the rela- 
tion between galaxy mass and metallicity (e.g., Lequeux et al. 
1979; Tremonti et al. 2004; Lee et al. 2006). Since flie metal 
content of NLRs is the result of the past star-formation his- 
tory in the host galaxies, the metallicity of galaxies and that 
of NLRs should be closely related. Galaxies are characterized 
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Fig. 4. The flux ratios of HzRGs plotted on the diagnostic diagram on C iv/He ii versus C iii]/C iv. Our new data are shown with 
red-filled circles, while the data compiled by Nagao et al. (2006a) are shown with blue-open circles and arrows. The green-filled 
star shows the line ratios measured in the composite spectrum of the 9 HzRGs observed by us. The shock models by Allen et al. 
(2008) are also plotted for comparison. Gray-solid lines denote the predictions of pure shock models with magnetic parameters of 
Bln^l^ - 0, 1, 2, and 4 yuG cm-'^^, solar abundance, and shock velocities in the range 200 km s"' < v, < 900 km s"' (constant 
velocity sequences of 200 km s ' and 900 km s ' are denoted by gray-dashed lines). The gray-dotted line denote the prediction 
of shock plus precursor models with Z?/n'^^ = 2 //G cm^**^, solar abundance, and 200 km s ' < Vj < 400 km s Black-solid and 
black-dashed lines are the same as those in Figure 3. 



by a well defined mass-metallicity relation even at high red- 
shift (Maiolino et al. 2008). On the other hand, AGN luminosi- 
ties and host galaxy masses should also be closely related, if 
the Eddington ratio is roughly the same within this class of ob- 
jects and the correlation between the mass of supermassive black 
holes (SMBHs) and that of galaxy spheroidal components holds 
also at high redshift. Therefore, by assuming a narrow range of 
the Eddington ratios and the MeH-galaxy mass relation, the pos- 
itive correlation between NLR metallicity and AGN luminosity 
is naturally expected. Indeed, some observational studies report 
that the Eddington ratio of high-z quasars is limited in a narrow 
range (e.g., Kollmeier et al. 2006; Trump et al. 2009). 

However, this scenario has a serious problem. We found 
that the NLR metallicity shows no significant redshift evolution, 
up to z ~ 4. On the contrary, the mass-metallicity relation in 
galaxies shows a significant redshift evolution, at least in the 
< z < 3 range (e.g., SavagHo et al. 2005; Erb et al. 2006; 
Liu et al. 2008; Maiolino et al. 2008). If the relation between 
NLR metallicity and AGN luminosity is simply caused by the 
galaxy mass-metallicity relation, the NLR metallicity is also ex- 
pected to show a significant redshift evolution, in contrast with 
our observational results. Below we discuss possible reasons for 
this apparent contradiction and the implications of the relation 
between NLR metallicity and AGN luminosity. 



A possibility to reconcile this problem is that the apparent 
lack of evolution is the extreme consequence of galaxy downsiz- 
ing on metal enrichment. We know that low mass galaxies evolve 
slowly, even from the chemical point of view, on a prolonged 
time scale extending to the current epoch. In contrast, massive 
galaxies reach their chemical maturity on short time scales and 
at high redshift (Maiolino et al. 2008 and references therein). 
Since the host galaxies of HzRGs are very massive (e.g. De 
Breuck et al. 2002), their chemical evolution may be completed 
at much earlier epochs than observed in the current sample, i.e. at 
z > 3. This scenario predicts that the NLR metallicities in lower- 
luminosity HzRGs (hence probably hosted in less massive hosts) 
may show evolution even at z < 3. This prediction can be tested 
with sensitive spectroscopic observations of faint HzRGs. 

Alternatively, the relation between NLR metallicity and 
AGN luminosity may be independent of the mass-metallicity re- 
lation of galaxies, since the Eddington ratio is likely not univer- 
sal. Recently it was reported that the BLR metallicity in quasars 
is correlated with the Eddington ratio (e.g., Shemmer et al. 2004; 
see also Nagao et al. 2002b; Shemmer & Netzer 2002). This sug- 
gests that the NLR metallicity may be also correlated with the 
Eddington ratio, although the physical origin is not clear In this 
case, the correlation between the AGN luminosity and the NLR 
metallicity is independent of the galaxy mass-metallcity relation. 
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Table 3. Averaged Diagnostic Flux Ratios of HzRGs 



Sample Number <z> < logL(Heii) > Znlr/Zq Civ/Heii Cra]/Civ 



1.2 < z < 2.0 


41.5 < logLCHeii)* 


< 42.5 


8 


1.70 


41.96 


0.80+"8 

-0.50 


1.82+!!-^' 

-0.37 


0.77+0;' 

-0.16 




42.5 < logLCHeii)* 


<43.0 


7 


1.66 


42.70 


1 97+0.34 
^^•^'-0.33 


1 79+0.41 
^ - '-^-0.33 


49+'"' 




43.0<logL(Heii)* 


<45.0 


2 


1.81 


43.32 


1 36+°'^' 


98+''-''* 

"•^"-0.27 


n QQ+o 40 

"•'^^-0.28 


2.0 < z < 2.5 


41.5 < logL(Heii)* 


< 42.5 


4 


2.28 


42.25 


77+0''' 

' -0.77 


1 34+°-" 

^•-^ -0.40 






42.5 < logL(Heii)* 


<43.0 


10 


2.28 


42.78 


1 54+''-3'* 

^•-^^-0.31 




65+° " 




43.0 < logLCHeii)* 


<45.0 


6 


2.24 


43.22 


1 53+''-2s 


1 QO+0.22 


56+" 


2.5<z<3.8 


41.5 < logLCHeii)* 


< 42.5 


4 


2.78 


42.39 




1 47+0.24 


94+""' 




42.5 < logLCHeii)* 


<43.0 


9 


3.06 


42.73 


1 ^ 9 +0.25 
'-•■-^ -0.22 


1 46+''-2' 


61+"'" 




43.0 < logLCHeii)* 


<45.0 


7 


2.67 


43.40 


1 64+''-23 

^•"^-0.24 


1 30+""' 


52+" "' 

"•-^ -0.08 



" Assuming the hydrogen density of iin = lO'^cm 
* Line luminosity of He ii/ll640 in units of ergs s"' . 
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Fig. 5. Flux ratios averaged logarithmically for each redshift 
bin, compared with the model predictions on the C iv/He ii ver- 
sus C iii]/C IV diagram. The dotted bars denote the RMS of the 
data distribution, and the solid bars denote the estimated errors 
on the averaged values. 



Whatever the scenario, our study strongly suggests that 
HzRGs completed their major chemical evolution in the very 
high-z universe, i.e., z > 4. If the minimum time-scale for a sig- 
nificant enrichment of carbon (~ 0.5 Gyr; e.g., Matteucci 2008) 
is taken into account, the major epoch of the star formation in 
HzRGs may have occurred at z > 5. 

Note that our conclusions are in contrast with some of the 
earlier studies mentioned in section 1 . The most significant dif- 
ference between previous works and our own is the adopted di- 
agnostics, i.e., we do not use N v emission to investigate the NLR 
metallicity. While the N v emission is too weak in most HzRGs 
to be measured accurately, all of the lines used by us (Civ, 
Hen, and Cm]) are so strong that we can measure their fluxes 
rather easily. Therefore we can investigate the chemical evolu- 
tion of HzRGs without relying on upper-limit data. Another se- 
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Fig. 6. Flux ratios averaged logarithmically for each luminosity 
bin, compared with the model predictions on the C iv/He ii versus 
C iii]/C IV diagram. The error bars are the same as in Figure 5. 



rious problem in using the Nv emission in metallicity studies 
is that the flux ratios including N v are sensitive not only to the 
NLR metallicity but also to the ionization parameter as discussed 
in the following section. 

5.3. The Nv emission in HzRGs 

Given the high S/N of our spectra, the N v emission line is se- 
curely detected in 5 out of 9 objects. Nitrogen is a secondary 
element and thus its abundance is expected to be proportional to 
the square of the metallicity [i.e., N/O oc O/H, or equivalently, 
N/H oc (O/H)^]. Therefore, the Nv line is widely regarded as a 
good metallicity indicator, at least for BLRs (e.g., Hamann & 
Ferland 1992, 1993; Hamann et al. 2002; Dietrich et al. 2002, 
2003; Nagao et al. 2006c). Accordingly the N v line is sometimes 
regarded as a metallicity indicator also for NLRs (e.g., van Ojik 
et al. 1994; De Breuck et al. 2000; Vernet et al. 2001 ; Overzier et 
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Redshift z 

Fig. 7. The NLR metallicity versus redshift diagram, plotting 
metallicities of three subsamples with log L(He ii) = 42.5 - 43.0 
at redshift range as follows: 1.2 < z < 2.0, 2.0 < z < 2.5, and 
2.5 < z < 3.8 (see also Table 3). These metallicities and the er- 
rors are estimated from the C iv/He ii versus C iii]/C iv diagram. 
Transverse bars show redshift ranges. 




43 44 

logZ(Hell) 

Fig. 8. The NLR metallicity versus luminosity diagram, plot- 
ting metallicities of three subsamples at luminosity range as fol- 
lows: 41.5 < logL(Heii) < 42.5, 42.5 < logL(Hen) < 43.0 
and, 43.0 < logL(Hen) < 45.0. These metallicities and the er- 
rors are estimated from the C iv/He ii versus C iii]/C iv diagram. 
Transverse bars show luminosity ranges. 



al. 2001; Humphrey et al. 2008). We thus examine the Nv/Civ 
and N v/He ii flux ratios as a function of redshift, to infer pos- 
sible redshift evolution of the NLR metallicity independently of 
the diagnostic diagram of C iv/He ii versus C m]/C iv. These two 
flux ratios are plotted as a function of redshift in Figure 9. Our 
new data at z > 2.7 clearly show the existence of HzRGs with 
high Nv/Civ and Nv/Heii flux ratios even at z > 2.7, which 
seems consistent to the idea that there is no significant evolution 
in the NLR metallicity of HzRGs as suggested by the C iv/He ii 
and C iii]/C iv flux ratios. 

Figure 9 shows, however, the presence of a large number 
of HzRGs with a very weak N v emission. In particular, some 
HzRGs show N v/He II < 0.1 and Nv/Civ < 0.1 while others 
show N v/He ii ~ 1 and N v/C iv ~ 1 . With the models presented 
by Vernet et al. (2001), these flux ratios would correspond to 
Znlr < OAZq for Nv-weak HzRGs and Znlr ~ 4Z0 for Nv- 
strong HzRGs. This would suggest that the NLR metallicity of 
HzRGs varies from object to object in a very wide range (> 1 
dex). This seems inconsistent with the results obtained through 
the analysis of the C iv/He ii and C iii]/C iv emission-line ratios 
(Figure 3). To investigate the origin of this apparent contradic- 
tion, in Figure 10 we plot the observational data on the N v/He ii 
versus N v/C iv diagram and compare them with our model cal- 
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Fig. 9. The flux ratios of new sample plotted on diagram of 
N v/He II (upper) and Nv/Civ (lower). Our new data are shown 
with the red-filled circles, and the previous data compiled by 
Nagao et al. (2006a) are shown with the blue-open cuxles and 
arrows. 



culations. The distribution of the observational data on this di- 
agram seems consistent with the area spanned by model grids. 
The near horizontal lines of constant- f/ in Figure 10 suggest 
that the N v/He ii flux ratio is sensitive to the ionization parame- 
ter but insensitive to the NLR metallicity. The N v/C iv ratio ap- 
pears to depend on both U and Znlr- In other words, the increase 
of the ionization parameter for a given metallicity results in high 
N v/He II and N v/C iv flux ratios, while the increase of the metal- 
licity for a given ionization parameter results in high N v/C iv 
but nearly constant Nv/Heii flux ratios. This suggests that the 
N v/He II flux ratio is sensitive to the ionization parameter rather 
than to the NLR metallicity. This is demonstrated more clearly 
in Figure 11, where we show the Nv/Heii versus Ciii]/Civ di- 
agram. The narrow coverage of the model grids in this diagram 
suggests that both N v/He ii and C ni]/C iv flux ratios depends 
mainly on the ionization parameter but are nearly independent of 
the NLR metallicity. Note that the observational data appears to 
be distributed in a wider range with respect to the one predicted 
by the models. Indeed, only the models with riu = 10"* cm"^ 
are plotted in the diagram and most of the observational data are 
upper limits. 

The physical reason why the N v/He n flux ratio is less sen- 
sitive than the N v/C iv flux ratio is shown in Figure 2 and de- 
scribed in section 4. Both N v and C iv are collisionally excited 
emission lines and thus become weak at high metallicity where 
the equilibrium temperature is low. However, N v emission de- 
creases much more slowly with metallicity than C iv, because 
nitrogen is a secondary element. This is the main reason of the 
metallicity dependence of the N v/C iv flux ratio. However, es- 
pecially at U > 10"^ *', the metallicity dependence of the Nv 
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Fig. 10. The flux ratios of HzRGs plotted on the diagnostic di- 
agram on N v/He ii versus N v/C iv. Red-filled circles, blue-open 
circles, arrows and green-filled star are the same as those in 
Figure 3. Solid and dashed lines are the same as those in Figure 
2. 




CIVX1549/HeIIX1640 



Fig. 12. The flux ratios of HzRGs plotted on the diagnostic dia- 
gram on N v/C IV versus C iv/He ii. Red-filled circles, blue-open 
circles, arrows and green-filled star are the same as those in 
Figure 3. Solid and dashed lines are the same as those in Figure 
2. 
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Fig. 11. The flux ratios of HzRGs plotted on the diagnostic dia- 
gram on N v/He n versus C iii]/C iv. Red-filled circles, blue-open 
circles, arrows and green-filled star are the same as those in 
Figure 3. Solid and dashed fines are the same as those in Figure 
2. 



emission is small (see the top panel of Figure 2). As already 
mentioned in section 4, the He ii emission is also insensitive to 
the metallicity because it is a recombination line. This is demon- 
strated in Figure 12, where we plot the N v/C iv versus C iv/He ii 
diagram. Here ZNLR-sequences are mostly perpendicular to U- 
sequences, and ZNLR-sequences go towards the upper-left direc- 
tion with increasing NLR metallicity. However this diagram is 
less useful than the C iv/He n versus C iii]/C iv diagram to infer 
the NLR metallicity, because N v/C iv flux ratios of HzRGs are 



difficult to be measured accurately; indeed most of the observa- 
tional data in Figure 12 are upper-limit values. 

In Figure 9-12, the inferred ionization parameter for the 
N v-detected HzRGs is much higher than that suggested by the 
diagram in Figure 4. This is due to the higher ionization poten- 
tial of the N'*+ ion (77.4 eV) compared to the C^+ (47.9 eV) and 
He^^ (54.4 e V) ions; i.e., the N v emission tends to arise from re- 
gions characterized by higher ionization parameters with respect 
to C IV and He ii. Although multi-zone photoionization models 
for NLRs (e.g., Ferguson et al. 1997a; Nagao et al. 2001b) solve 
this discrepancy in ionization parameter, it is beyond the scope 
of this paper and thus they are not discussed further. 



6. Conclusion 

To investigate the metallicity evolution of galaxies in the high-z 
universe, we have performed new spectroscopic observations of 
9 HzRGs at z > 2.7, which complement data available at lower 
redshifts. By comparing the total set of data of 57 HzRGs at 
1 ^ z ^ 4 with photoionization models, we found the following 
results. 



Our analysis of the emission-line flux ratios of C iv. He ii and 
C III] suggests that there is no significant chemical evolution 
in the redshift range of 1 ^ z 5 4. 

We found a positive correlation between the NLR metallicity 
and the AGN luminosity. There are two possible origins for 
this correlation: the galaxy mass-metallicity relation and a 
dependence of the NLR metallicity on the Eddington ratio. 
The non-evolution of the gas metallicity of the NLRs implies 
that the major epoch of star formation in the host galaxies is 
atz > 5. 

We detected the N v emission in the spectra of 5 HzRGs at 
z S 2.7, and found that there are some HzRGs with high 
N v/C IV and N v/He ii flux ratios even at such high redshift. 
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5 . However, high N v/He ii ratios do not always mean high NLR 
metallicities, but correspond mainly to high ionization pa- 
rameters. 
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